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r~^ ■ Abstract: We calculate the scale dependence of the bispectrum and trispectrum m (quasi) local 

t~^ , models of non-Gaussian primordial density perturbations, and characterize this scale dependence in 

^^ I terms of new observable parameters. They can help to discriminate between models of inflation, 

since they are sensitive to properties of the inflationary physics that are not probed by the standard 
observables. We find consistency relations between these parameters in certain classes of models. 



o 

f~^ ' We apply our results to a scenario of modulated reheating, showing that the scale dependence of 

non-Gaussianity can be significant. We also discuss the scale dependence of the bispectrum and 
trispectrum, in cases where one varies the shape as well as the overall scale of the figure under 
(^j consideration. We conclude providing a formulation of the curvature perturbation in real space, which 

%_] . generalises the standard local form by dropping the assumption that /nl and (Jnl are constants. 
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1. Introduction 

Inflation is the simplest framework which explains the origin of the observed power spectrum of 
temperature fluctuations in the cosmic microwave background |l|. It is now widely accepted that 
non-Gaussianity is a powerful probe to discriminate between the many currently viable inflationary 
models ||, ||, ^, [s], ^, |^, || . It is usually parameterized in terms of a single constant parameter, /nl, 
corresponding to the amplitude of the bispectrum normalized to the square of the power spectrum of 
primordial curvature fluctuations M lOl . More recently it has become common to further characterize 
local non-Gaussianity including the two non-linearity parameters associated to the trispectrum, called 
5nl and tnl, again treating them as constants |ll|. However, it has been recently pointed out, both 
from theoretical ||lj, ^ |lj, ^ and observational viewpoints [Q, that /nl is not necessarily constant. 
We show the same holds true for (^nl and tnl ■ As happens with the power spectrum and the spectral 
index, they are characterized by a scale dependence, that which denote respectively with 7^/NL , "-^nl 
and n^-j^L ■ For example, if /nl is large and positive on large scale structure scales [^ ^ ^ ^ , but 
has a smaller value on the largest CMB scales then this would require that /nl is scale dependent. Any 



scale dependence of the non-linearity parameters provides a new and potentially powerful observational 
probe of inflationary physics. 

In this work, we discuss a new approach to study the scale dependence of the non-linearity pa- 
rameters, based on the JiV- formalism pi], p2l p3|. This allows us to obtain an expression for the 



curvature perturbation that generalizes the local Ansatz introduced in M, 10 , and that contains the 
aforementioned scale-dependent parameters. For the single field case, the curvature perturbation can 
be schematically written as 

Ck = C^ + ^/Sl (1 + "/NL In k) (C° * C°)k + ^5&L (1 + 'W In k) (C° * C'' * C'')k , 

where (^ is a Gaussian variable, and /^^ ^nd g^^ ^re constants. Our approach allows us to directly 
calculate j^/^l j '^-^nl ^-nd ^^tnl in models with an arbitrary inflationary potential and an arbitrary 
number of fields, assuming slow-roll inflation. We also assume the field perturbations are Gaussian at 
Hubble exit. Our results depend only on the slow-roll parameters evaluated at Hubble exit, and on 
the derivatives of N^ the e-folding number. In particular, we find that n/^^ ^nd rig-^^ are sensitive to 
third and fourth derivatives of the potential along the directions in field space that are responsible for 
generating non-Gaussianities. These do not in general coincide with the adiabatic direction (during 
inflation) and such features cannot therefore be probed by only studying the spectral index and its 
running p3] . In the case that a single field generates the curvature perturbation there is a consistency 
relation between ri/^L and jt-tnl which is the derivative of the consistency relation between /nl and 
tnl- We explicitly show how this consistency relation is violated in multiple field models. 

In the framework of slow-roll inflation, there are various ways to generate large non-Gaussianity, 
in models in which more than one field play a role during the inflationary process. This is the case of 
multiple field inflation H |6[ |27[ |28[ |2|, ^, |l|, H g |4| , in which two or more fields contribute to 
the curvature perturbations. But there are also approaches in which, although more than one field is 
light during inflation, only one of them contributes significantly to the curvature perturbations (the 
most studied examples are the curvaton [Q and modulated reheating [^6[ ^ scenarios). In this 
work, we apply our general findings to this last class of models. We consider set-ups in which an 
isocurvature field remains subdominant during infiation (as required in order to have an observable 
level of non-Gaussianity pq]), but represents the main source of curvature fluctuations after inflation 
ends. In this case, neither the spectral index of the power spectrum of curvature perturbations, nor 
its running are sensitive to the third and fourth derivative of the potential of the subdominant field. 
Hence the scale dependence of non-Gaussian parameters provide a unique opportunity to probe self 
interactions in these scenarios. As an example, in the modulated reheating scenario, it is possible 
for any of the non-Gaussian parameters to be large. We show that if the modulaton field has self 
interactions, for example a quartic potential, then all of /nl, "-/nl, 5nl and rig^j^ can be large and 
provide novel information about the mechanism which generates curvature perturbations. We will also 
consider mixed scenarios in which the infiaton perturbations are not neglected pfl, Eol Efl, E2] . We note 
that the scale dependence of equilateral type non-Gaussianity is also of theoretical and observational 
interest ||, ||, ||, |6 . 



We have previously shown ||1J] that provided one scales all three sides of the bispectrum at the 
same rate then tt-z^l is a constant (and hence it is simplest to focus on an equilateral configuration). 
We show a similar result for the trispcctrum parameters. Since it may be of interest to consider 
more general variations in which one changes the shape of the figure under consideration, we also 
consider this case. We find the combination of shape and scale dependence which maximizes n/^^ 
and show that it is never significantly larger than the standard result, in which one keeps the shape 



fixed. However we single out interesting limits in which there is no scale dependence, corresponding 
to squeezed figures. 

While in most of the paper we work in momentum space, in the last part we also discuss how to 
describe our results in coordinate space. We provide an expression for the curvature perturbations 
in real space, that generalizes the simplest local Ansatz [|0| , and that exhibits directly in coordinate 
space the effect of scale dependence of non-Gaussian parameters. 

The plan of our paper is as follows: In Sec. y we extend and simplify the results from our previous 
paper to give general results for the non-linearity parameters, including those which measure the 
trispectrum. In Sec. || we reduce the results to general single field models and derive a consistency 
relation. In Sec. we consider simple one or two-field models in which one field, e.g. the curvaton, 
generates non-Gaussianity but we do not exclude the Gaussian perturbations from the infiaton. Many 
popular models in the literature fall into this class and the reader may choose to skip straight to 
this section where the results and notation are significantly simpler. As an explicit example we study 
modulated reheating. In Sec. || we consider in detail the scale dependence of the bispectrum and the 
trispectrum, and how this can be affected by the shape of the triangle or quadrilaterum. In Sec. O we 
consider how to generalize the coordinate space expression of the curvature perturbation to include 
scale dependence. Finally we conclude in Sec. m 

2. General results 

In this section, we discuss a new approach to analyze the scale dependence of quasi-local non- 
Gaussianity, by means of a suitable implementation of the 5N formalism. Using the 6N formalism 



21, ^ ^, the curvature perturbation for a system of n scalar fields (jia is given by the expression 

C(i/,x) = J2Naitf,U)Sr{U,^) + i^iV,b(i/,i0^r(i.,x)<5/(i„x) + • • • , (2.1) 

a ab 

where tf labels a uniform energy density hypersurface and U denotes a spatially flat hypersurface. The 
result is valid on super-horizon scales where spatial gradients can be neglected. In this work, we do 
not consider secondary effects on curvature perturbations arising from late-time physics (see P7| for 
a review), nor the effects of possible isocurvature modes during the late universe. The quantities Na 
and Nab denote derivatives of the number of e-foldings along the scalar fields. We choose U as a time 
soon after the horizon crossing of all the modes of interest. Written in momentum space, Eq. ( p.l| ) 
reads 

Ck(i/) = 5]A^a(i/,i.)^0£(iO + i^iVa6(i/,i.)(<50"(iO*'5/(iO)^ + --- . (2.2) 

a ab 

Here k < a{ti)H{ti), since we focus on super-horizon scales, and 7k- denotes a convolution: 

(srik) * scb\u))^ ^Jj0f mtos^i-<,(t^) ■ (2-3) 

To analyze the statistical properties of the curvature perturbation it is useful to express the results 
in terms of scalar perturbations evaluated at horizon crossing S(f>^{tk). Assuming the fields (j)a obey 
slow roll dynamics during inflation and have canonical kinetic terms, S(J3^(tk) are Gaussian up to slow 
roll corrections p^, E^. These corrections are irrelevant in cases where the non-Gaussianities are 
large, |/nl| >> 1 or |.gNL| ^ 1- Therefore, in our analysis we take the fields S(j)^{tk) to be Gaussian 
at horizon crossing, k = a{tk)H{tk)- In Q the result (2^) was expressed in terms of S(j}l.(tk) by 



setting ti — )■ tk{k). This makes the coefficients Nat... imphcitly dependent on k. In this work, we 
follow a different approach, choosing a fixed U for all observable k modes, and explicitly solving for 
the perturbations at U as a function of S(f>^(tk). Besides being more transparent, this method allows 
us to easily write down explicit results for the scale-dependence of non-linearity parameters. The two 
approaches are compared in detail in Appendix pi 

We first note that, assuming slow roll, the evolution of super- horizon scale fiuctuations 6(j)"'{t,x.) 
from some initial spatially flat hypersurface at to < ti to the spatially fiat hypersurface at ti can be 
expressed in terms of the Taylor expansion 

Here we have also assumed the fields have canonical kinetic terms, i.e. the metric in field space is fiat. 



The result (2.4) follows directly from the application of the 6N formalism where any super-horizon 
region, labeled by x, evolves as a separate FRW universe with its own initial conditions. Since we 
assume that slow roll conditions are satisfied, the initial conditions are set by the field values {'/>'' (^o)} 
alone, i.e. any dependence on the field time derivatives can be neglected. Therefore, 

sr{u,^)^riu){{^''{to) + 6c^''{to,^)})-r{u){w\to)}) , (2.5) 

where <j)°-(ti){...) denote FRW solutions with the initial conditions set at ip- Eq. ( p!4| ) is obtained 
by expanding this with respect to {(50''(io,x)} while keeping fixed the number of e-foldings between 
to and ti. This corresponds to choosing to and ti as spatially fiat hyper-surfaces, since it amounts 
to comparing different realizations of FRW universes that all undergo the same number of e-foldings 
between to and ti. 



The coefficients in Eq. (2.4) can easily be computed by solving the slow roll equations of motion, 
iH(j)a = -Va. We find 

riU) = r{ta) ~ V2^Ha^/ao) + O (^e^/^\n^{a,/ao)^ , (2.6) 

where the slow roll parameters are evaluated at to and defined as usual: Eq == (T4/(3-ff^))^/2 and 
Vab = Vab/iiH"^) (with A/p = 1). In the following we neglect the slow-roll suppressed corrections 
©(e'^/^ln (ai/ao)), where C'(e'^/^) denotes terms involving powers of Eq and i^ab up to 3/2. The validity 
of this approximation is discussed in more detail below and also in Appendix [Aj. Differentiating 
Eq. (^^) once with respect to the initial field values, and keeping \n{ai/ao) fixed, we find 

d(j)°-iti) 

cij,br.\ = ^ab + eab ln{ai/ao) , (2.7) 

09 {to) 

where we have defined 

(2.8) 



The higher order derivatives in Eq. (2.6) can be computed in a similar way and the results are given 
in Appendix ^. 



By substituting Eq. (2.4) into the coordinate space expression for the curvature perturbation ([2.1| 



taking the Fourier transform and thereafter setting to = t/,, we arrive at the result 

Ck - E Ck '" + E /«^(fc)(c'''" * c^''')k + E 5abc(fc)(c^''^ * c""'' * c^'^)k + • • • . (2.9) 

a ab abc 



aG,o 



Here Ck ' are Gaussian fields defined as 



Sab + £ab In 



fli-ffj 



{aiHi - k) 



(2.10) 



The Gaussianity of this quantity follows from our assumption of the perturbations 5(j)^{tk) being 
Gaussian at the horizon crossing. For brevity, from now on we suppress the time arguments of the 
derivatives of iV, denoting Nab... = Nab...{tf,ti). The theta function in Eq. (2.1C) is included to 
constrain the convolutions to only include super- horizon scales, k < aiHi. The matrices fab{k) and 
9abc{k) are given by 



Iab{k) 

gabcik) 



1 Nab 



^E 



2 NaNb 2 

1 Nabc 

6 NaNbNc GNaNbNc 



NaNb k 

1 



E 

d 



(SNda 



F, 



(2) 



dbc 



^.^il)ln^ 



(2.11) 
(2.12) 



where k < QjHj and F, 



(m) 



abi 



denotes the fc-independent part of the m:th order coefficient in Eq. (p.4| 



They are proportional to combinations of slow roll parameters and their explicit expressions are given 
in Appendix |A|. 

Our results are derived to first order in InaiHi/k. In Eqs. ( 2.11 ) and ( 2.12| ) the terms proportional 
to InaiHi/k represent small corrections to the fc-independent parts except in the cases where fab{ki) 
and gabc{ki) are comparable to slow roll parameters or even smaller. Such components, however, do 
not generate observable non-Gaussianity and therefore do not play an important role in our discussion. 
For this reason, we can safely perform the expansion in hiaiHi/k. Since the higher order terms arising 
in this expansion are further slow roll suppressed, and since we can choose ti such that the logarithms 
never get larger than 0(10) for the super-horizon modes in our observable universe, we can truncate 
the expansion at first order. 

Instead of expanding in \naiHi/k, we can also choose one of our observable super-horizon modes 



as a pivot-scale, fcp < aiHi, and expand Eqs. (2.11) and (2.12) around this point. To first order in 
ln(fc/fcp) the results are given by 

(2.13) 

(2.14) 

(2.15) 



fab{k) = fab{kp) 1 -I- n/,afc In— , 



Jabc{k) = gabc{kp) 1 +ng,afcc In^ 



where we have defined^ 
_ d\n\fab\ _ 



nf^ab 



dink 



^ Nab 



^g,abc 



dln\gabc\ 
dink 



Nda^ p(2) ^ Nd p(3) 
Nabc 



Z^[^ AT . ^dbc 



Na. 



be 



dabc 



Provided that nf^ab and rig.abc are not much larger than 0(0.01), truncating the above series at first 
order leads to an error of a few per cents at most. Neglecting slow roll corrections, we can write 
fab{kp) = Nab/C^NaNb) and gabc{kp) = Nabc/{QNaNbNc). This precision suffices when treating the 
fc-independent terms in our expressions, since we are only interested in computing scale-dependencies 



^It is important to realise tliat n-f^b ("-g.abc) is only defined in the case where f^b J^ {9abc ¥" 0) in th^ limit 
k —> aiHi. If fab (dab) vanishes, it is convenient to define the derivative in Eq. (2.15) to be identically zero. 



to leading order in slow roll. In what follows, we will therefore always write the constant terms to 
leading order precision in slow roll. 

Finally, using Eqs. ( |2.13| ) and ( [j.l4| ), we can express Eq. ( ^.9[ ) as 



ab 



k 



+ E 5afcc(fcp) ( 1 + Ug^abc \n^ ) (C""'" * C"^'" * C°'')k 



abc 



This result is the starting point for our analysis of the scale-dependence of non-linearity parameters. 
Explicit expressions for njab and Ug abc are given in Appendix W and the scale-dependency arising 
from the fields C^. '° can be computed using standard methods. Therefore, using Eq. (2.16|) we can 



explicitly compute the scale-dependencies of /nl, .9nl and tnl for any model with slow roll dynamics 
during inflation. 

2.1 Two point function and power spectrum 

Here we re-derive some well-known results for the scale dependence of the spectrum of curvature 
perturbations; they will be useful in what follows for analyzing the scale-dependence of bispectrum 
and trispectrum. The two point function of the scalar field perturbations S(t>l^{tk) at horizon crossing 
is given by 

{SK,itk)S(bi^{tk)) ^ (2^)3j(ki +k2)|^ ( ^^^2^ '^^ ) [Sab + 2c{l-Sab)eab] , (2.17) 

where c = 2 — ln2 — 7~0.73 with 7 being the Euler-Mascheroni constant. Both the diagonal a — b 
and off-diagonal a ^ b components are given to leading order in slow roll. Note that although the 
off-diagonal components are slow roll suppressed compared to the diagonal components, their scale 
dependence has no further suppression and therefore gives a contribution comparable to the scale- 
dependence of the diagonal components. Therefore, we need to retain the off-diagonal contributions 
in our analysis. 



Using this together with Eqs. (2.10) and (2.16), we can express the power spectrum of C, defined 
by (CkiCk.) = (2^)3,5(ki +k2)P(fci), in the form 

0^2 

(2.18) 

K~ ft;- 

ab 

Here 



Vabik) 



P{k) = ^Vik) - 


ab 


^)^-- 


Sab 1 - 

. V 


k \ ( k 
- 2eH In 1— ) + 2ea6 c - In -— 

kp) V K 



(2.19) 



and we have defined eh = —H/H^ and c = c + ln(aj_ffi/A:p). Subleading slow roll corrections are again 
neglected in the scale-independent terms. 
Defining a quantity 

nab-l= ^^^ = ~6ab{2eH + 2eab) - l{l - S^b) , (2.20) 

we can write the spectral index as 



ab 

This agrees with the result given in ]2 



""i^- dink - ^\V J '""" ^ " "^" " EcN^ 



2.2 Three point function and /] 



NL 



We now proceed to apply our approach to derive the scale dependence of non-linearity parameters. 
Using previous definitions we can write /nl in a general multiple field case as 



/NL(fcl,fc2,fc3) 



B{ki,k2,k3) 



6P{ki)P{k2) + 2 perms 

5j2abcdikik2)~^Vaciki)rbd{k2)fcdik3) + 2pevins 
3 {kik2)-^r(ki)r{k2) + 2peTms 



(2.22) 



where the bispectrum is defined by (27r)'^(5(ki + k2 + k3)_B(fci, fc2, fca) — (CkiCkaCkj 
In the equilateral case, ki = k for i = 1, 2 and 3, this simplifies to 

^i:abcd'Pac{k)VM{k)fcd{k) 



hUk) 



V{kY 



and using Eqs. ( 2.13 ) and ( 2.19 ), we find the scale dependence of fni^ik) is given by 

_dln|/NL(fc)| 1 



din /c 



/NL(fcp) 



Here we have defined 



and 



,ab ^ 5 NgNbNab 



/ ^ /nl (2'^multi,a + ^ f ,ab) 



/NL(fcp)=^/5i 
ah 



N, 



'^multi.a = naa ~ "C ^ ^^(1 - 5ac)f-ac-rf 



Na 






cd 



EbN^ 



Nn 



(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2.27) 



All the quantities in Eq. (2.24) depend on combinations of slow roll parameters and on the constant 



coefficients Na,Nab in the 6N expansion (recall that the explicit expression for ri/,a6 is given by 



Eq. (A.7)). For a given model these can all be regarded as known quantities and the scale-dependence 
of /nl can therefore be directly read off from the above result without doing any further computations. 



In Eq. (2.24 ) we can clearly identify two sources of scale dependence. The contribution proportional 
to Uf^ab follows from the non-linear evolution of perturbations outside the horizon ||lj]. The part 
proportional to rimuiti,a is associated with the scale dependence of factors of the form Vac/'P in equation 
( [2.23| ). It is present only in the multi- field case (indeed for a single field model this factor is equal to 
unity) and arises due to the presence of multiple unrelated Gaussian fields C*^'" in the expansion of C, 
in (^.16). This generates deviations from the local form and makes /nl scale-dependent even if the 



perturbations would evolve linearly outside the horizon. Indeed, by setting nj^ab = we recover the 
results of a multi-local case analyzed separately in [|lj] . 

As shown in [Q, the scale dependence of /nl is given by the same result ( 2.24 ) for the class 
of variations where the sides are scaled by the same constant factor, ki — > aki. For such shape- 
preserving variations where only the overall scale of the triangle is varied, the result does not depend 
on the triangle shape. This holds at leading order in slow roll. Generic variations changing both the 
scale and the shape of the triangle are considered in Sec. ||. 



2.3 Four point function, ^nl and tnl 

The connected part of the four point correlator of (^ can be written in the form 

4 



(CkXk.CkaCkJ = {27:fSiJ2^,) 



TNL(fci, fc2, h,h, ku) ( P(fci)P(fc2)-P(|ki + kal) + 11 perm) (2.28) 



54 
25 



5NL(fci,fc2,fc3,fc4) ( P{ki)P{k2)P{k3) + Sperm' 



where we have defined fcy = jk^ + k^l. The functions tnl and ^nl are given by 



TNhih, k2, ks, ki, kis) 



gNL(fcl,fc2,fc3,fc4) 



^ {klk2ki3)-''J2abcdef'Pac{kl)Vte{k2)Vdf{kl3)fcd{k3)fef{k4) + H pCrm 

(fcifc2fci3)-3 75(^1 )7'(fc2)^(fci3) + 11 perm 
25 {kik2k3y^ J2abcdef Pa.d{ki)Pbe{k2)Pcf{k3)gdef{k4:) + 3 pcrms 
T (fcife/ca)-^ V{ki)V{k2)V{k3) + 3 perm ' 



,(2.29) 
(2.30) 



In the case of a square, k = hi (notice that tnl, but not 5nl, is sensitive to the angles between 
the vectors and different equilateral figures in general yield different results), the above expressions 
reduce to 



rNL(fc) = 4 ^ 

abcdef 



ffNL(fc) 



25 



E 

abcdef 



Vac{k)Vbe{k)Vdf{V2k)f,d{k)Uf{k) 
V{kY 

Vad{k)Vbe{k)Vcf{k)gdef{k) 



V{kf 



(2.31) 
(2.32) 



The scale-dependence can be computed similarly to the analysis of the bispectrum above. Using 
Eqs. (|1|), (|I§ and ( |2T9| ), we find 

nr^^ = "'^''j^';^^^' = TFT E ^Sl"" [(2nmuiti,a " ("c - 1) - 2eH)^bc - 2ebc + ^nj^abhc] ,(2.33) 



d\nk 

dln|gNL(fc)| 
din fc 



j,bcd 



/, N /_^ 3nL (3?^nmlti,a + ng^abc) 



g^hikp) 



where e// = -H/H^, and ||, |n| 



abed ^ NgNatN^dNd 

„,, _ 25 JVgiVbJV.JVabc 



Tnl 



e^p) = E ■^NL ''^fcc 



abed 
abe 



be 

5nl 



(2.34) 

(2.35) 
(2.36) 



The scale-dependencies are fully determined by the constant coefficients Na^Nab^Nabc in the 5N 
expression and by combinations of slow-roll parameters, which enter the results through Eqs. (A.7) 
and ( A. 8 ). Although the expressions appear lengthy in their general form, considerable simplifications 
typically occur when considering specific models. We will discuss examples in Sections and 0. 

Similarly to fi/^L ■, we can again distinguish two physically different contributions in the expressions 
for Ht-nl and %nl- The parts proportional to Uf^ab and Ug^abc in Eqs. (2.33) and (2.34) respectively 
arise from the non-linear evolution outside the horizon. The other contributions describe deviations 



from the local form due to the presence of multiple fields, similarly to what we discussed in the previous 
section. 

The results ( 2.33| ) and ( 2.34 ) hold not only for the special case of a square, but for any variations 
where all the sides are scaled by the same constant factor, k^ -^ aki. These variations preserve the 
shape of the momentum space figure and change only its overall scale. We will prove this result in Sec. 
g where we also discuss generic variations that simultaneously change both the scale and the shape. 

Having presenting our formalism and the general results, we will discuss in the next two sections 
applications to specific cases. 



3. General single field case 

We start by discussing models where the primordial curvature perturbation effectively arises from a 
single scalar field, which does not need to be the inflaton and we call cr. In this case, the functions 
/crcr(fc) and Qaaaik) appearing in the expansion of C, Eq. (2. IE), are up to numerical factors equal 
to f^hik) and gyii^{k), evaluated for the equilateral configurations. This can be seen directly from 
Eqs. (^.23|) and (^.32). We can therefore rewrite Eq. ( J.16|) as 



Ck = (^ 



^/NL(fc)(C° * C°)k + ^5NL(fc)(C° * C° * C°)k + 



(3.1) 



As we will discuss in Sec. ]^ this result applies for example to the curvaton scenario and modulated 
reheating in the limit where the inflaton perturbations are negligible. We therefore call all the models 



where the curvature perturbation can be expressed in the form (3.1) as general single field models. 

According to Eqs. (2.13) and (ET4), the non-linearity parameters /nlC^) and (7NL(fc) are now 
given by 



/NL(fc) = 

gNL(fc) = 



5 N" 

25 N'" 
54 ]V^ 



-nf^ 



In- 



1 ^ 

^ffNL In— 

/Cr 



where the primes denote derivatives with respect to a and ti/^l = 

in the Appendix [AJ, we obtain 



rif. 



explicit expressions (A. 7) and 
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Ig.aacr- Usiug the 

(3.4) 
(3.5) 



The same results can of course be directly obtained from Eqs. ( [2.24 ) and ( ^.34 ). If a is an isocurvature 
field during inflation, eg- = in the above expressions. 

For the general single field case Eq. ( |2.3l| ) further yields 



TNL(fc) 



/ 6/NL(fc) 



(3.6) 



up to scale-independent slow roll corrections. Therefore, the scale-dependencies of tnl and /nl are 
related by 

"tnl = 2n/„L ■ (3-7) 

This simple consistency relation is characteristic for general single field models. In multiple field 



models, the relation (5.6) is in general violated and consequently the result (3.7) is no longer true 



4. Two field models of inflation 

After considering single field models, in this section we discuss some scenarios in which more than 
one field can play an important role in the inflationary process. We focus on a class of models that 
contains the most important examples of inflationary set-ups characterized by large non-Gaussianity. 
Many models of inflation that generate sizeable non-Gaussianity are characterized by the presence 
of a field a, with significant non-Gaussian perturbations, that is isocurvature during inflation. The 
inflaton fleld also has its own perturbations, which for convenience can be considered as Gaussian. 
When the inflaton perturbations provide non-negligible contributions to the curvature fluctuation 
spectrum, the scenario is called a mixed scenario [^ ^ ^, ^. In order to generate large non- 
Gaussianity by means of the fleld a, it is required that a <^ (j), and hence Co- ^ £</, [Q. From this 
relation, it follows that the trajectory in fleld space while observable modes exit the horizon is nearly 
straight. Therefore it is a good approximation to treat the flelds as uncorrelated [BOl. We also make 
the common assumption that the potential is separable, 

W{a,(P)^U{(P) + V{(j) . (4.1) 

Hence, the only potentially non-negligible slow roll parameters in such a scenario are the following 

__H_ _U^ _V^ ^2 _ U'"U' ,2 _ y"'U' 

en - £0 — ^2 ' '^'t' ~ 3^2 ' ^^ ^ ^f£2 ' ^0 ~ QJJ4 ' ^<^ ~ QffA ■ ^^■^) 

In this case, the curvature perturbation reads^ 

C(k) = C^-^ + C^'^ + Uk) (C^-- * C^'-)^ + 5.(fc)(C^''^ * C'^ * en^^ ■ (4.3) 

Although the assumed form of Q is simplifled, in practice the vast majority of models in the literature, 
characterized by large quasi-local non-Gaussianity, satisfy the above Ansatz to a good enough accuracy 
for observational purposes. For this reason we will limit our attention to models with curvature 



perturbation satisfying Eq. (4.3) in this section. 

In the limit that f^ and g^ are independent of fc, we recover the multivariate local model |14|| . In 
the case that C,'^''^ = we have the general single fleld model we have analyzed in section ^ but here 
we assume this field was an isocurvature mode during horizon crossing. We will consider these two 
cases in more detail later in this section. 

The power spectrum is given by 

V^ik) = V^{k) + VAk) = V4k){l - w^k))-^ , (4.4) 

where we have introduced the ratio 

w^k) = ^ . (4.5) 

Note that neglecting all the slow-roll corrections, and hence also the scale dependence, w^ = N^/{N^ + 
N^)- To lowest order in slow roll, the spectral index n^ — 1 and tensor-to-scalar ratio tt satisfy the 



^We have used a simplified notation for this section compared to the rest of the paper. Since all cross terms such 
as P^cT ^re negligibly small in this scenario we use only a single index <j> or a where appropriate, e.g. for ri„ = r)„„ and 

5o" ^ QcJ(T(J- 
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following relations pll 



n^ - 1 = (ncr - l)wa- + (nil, - 1)(1 - w^) 

= -(6 - 4w^)eH + 2(1 - w„)t]^ + 2w„r^„ , (4.6) 

= 'El - 8 

= 8N-^l-w,), (4.7) 



where Vt ~ 8iJ^/(47r^) is the power spectrum of tensor perturbations and we have defined 

dlnVa , dlnV^ , dlnVr 



(4.8) 



The non-Gaussianity parameter /nl in the equilateral limit, and the trispectrum non-linearity 
parameters in the case of a square configuration, are given by 

/NL(fc) = ^wlik)Uk) , (4.9) 

rNL(fc) = 4wAk)^w,{V2k)f^{k) , (4.f0) 

25 
~9 



25 
ffNL(fc) = -u;:^(fc)g.(fc) . (4.ff) 



Therefore their scale dependence is given by 

_ dln|/NL| „, ^ , rfln|/a| ,. ^„. 



d In /c d In fc 



N„ ( V" , /I ^ ^/^^ 



4(l-u>,)(26H + r7.-r;^) + ^(— 2-\/2i^ry, (— -1) ), (4.13) 



n... =3K-"C) + 2^^J^ (4-14) 

= 6(l-w;a)(2eH + r/.-r/^) + iP( ^772-^/2eH?7^( — -1) I, (4.15) 



2N^ ( V" , / 1 ^ ^'"^^ 



3ff2 



UJn 



dlnfc 



^SNL = 3(n^ - n<;) + (4.16) 



3 A^ y " 

6(l-^,)(26H + ?7.-'?0) + ^^^ (4.17) 



^U^-3,. + V2i^^ --1 



1/2N 



where we have used the results derived in Sec. and the fact that Nss, N^a and their derivatives 



4>{T 



are negligible in the class of models we are considering, see Eq. (4.3). The quantities on the right 
hand side of each equation should be evaluated at an initial time ti shortly after the horizon crossing 
time of all the modes of observational interest. Observe that our results for Uf^-^ and Ug^-^ depend 
on the derivatives of the potential, in combinations that do not correspond to traditional slow-roll 
parameters. This turns out be useful to probe these quantities, that cannot be tested by the power 
spectrum and its derivatives. We are going to discuss this in detail in what follows. 

Observational constraints on the bispectrum are given in |l| while constraints on (/nl are given in 



52, B3l (see also M) and for both ^nl and tnl in [p5[- Forecasts for future constraints on all three 
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parameters are given in |5^, ^]j while forecast constraints on ti/^l are given in |16[| . There are currently 
no forecasts for how well the scale dependence of the trispectrum parameters could be constrained 



or measured. Observational constraints on a model with the form (4.3), without considering the 
scale-dependence of /nl or w^ , are given in [p8| . 

4.1 Limiting cases 

After presenting the general formulae for the two-field case, we discuss important examples of general 
single field inflation, that arise as limiting cases of the previous discussion of two-field inflation. 

4.1.1 Isocurvature single field 

In the case that a single field a, which is subdominant to the inflaton during inflation, generates the 
primordial curvature perturbation, one has w^ = 1 which implies Ua^ = n, rx — and No- 3> 1. 

In this scenario, it is useful to express the spectral index and its running up to second order, to 
understand which parameters are currently constrained by observations. From pO| , we have 

22 ^ \ o Zq/^ .\ I 2 



nc - 1 = -2eH + 2rj,+ \-— + 8c j e^ + -rf^ + ( 3 " ^c 1 e^V^ + ( 3 " ^c 1 ewia , (4.18) 
"C = -84 + 4eHf70 + 4eH77<T , (4.19) 

where c = 2 — ln2 — 7~ 0.73. Notice that, in the previous formulae, the slow-roll parameters £,a 
and ^.ip do not appear in the running of the spectral index, because they are weighted by negligible 
quantities. This implies that third and higher derivatives of the potential do not enter in the previous 
quantities. 

The non-Gaussianity observables (which follow as special cases of the formulae discussed in the 
first part of this section, when taking the limit Wo- — >■ 1) are 

, 5 5iV,, 25 25 7V,,, 

/NL=3/.-g^, 3NL = y5.-^^^, (4.20) 

_ rir^, ^ iV, V" ^, 5 sgn(JV,) [7^ V" 



^ iV„ V" , iV, fV"" 

'^SNL — '^ AT OU2 + I\r OIT2 ""/cr 



3 5NL V 8 3i?2 54 5NL 8 3if2 g^^ f«^ 54 5nl 67r2 ^ ' 

Here /nl and (7nl denote the non-linearity parameters evaluated for equilateral configurations at some 
pivot scale, k = kp. As discussed in Sec. g, fcp can be chosen as any of the super- horizon modes in 
our observable universe and the results are independent on this choice, up to subleading slow-roll 
corrections. 

In the previous formulae, we have presented several different ways of expressing ?t./nl and tt-^nl (in 



Eq. (4.23) we have dropped the negligible contribution rj^rT < 10"^). This is in order to make it easier 
to estimate their magnitude in different ways, depending on the available quantities. We also note 
that in some cases the previous formulae might include terms at different orders in slow roll, in which 
case one should neglect the subleading terms (since additional terms at the same order might also have 
been neglected). In general, they are suppressed by some combination of the tensor-to-scalar ratio, 
divided by non-linearity parameters. But their size could be significant, if a has either a large cubic or 
quartic self interaction. As we mentioned earlier, the power spectrum does not contain information on 
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these parameters, even if the running of the spectral index can be measured. Hence h/^l appears to 
be the best way of probing the cubic self interaction, while in principle rig^^^ could probe the quartic 
derivative of the isocurvaton field. 

Although we have written Uf^^-^ ^ 1//nl, the prefactor to 1//nl will in general depend on some 
of the same model parameters as /nl so one should not view the two parameters as being inversely 
proportional (an explicit example is given in [12|). In the case that /nl follows an exact power law 
behavior, /nl oc /c"-''nl then /nl and n/j^^ are of course independent. However if /nl = Aln(fc) + B 
where A and B are constants, then n/„L = A/(^ln(fc) + B) = A//nl- In this case rif-^-^ and /nl are 
not independent. Nonetheless one can easily check that the running of n/„L satisfies a/^j^ = —nl so 
it is a good approximation to treat rif^^ as constant provided that I^/nlI ^ 1- 



Consider, as a first example, the curvaton scenario [35 in the pure curvaton limit. In this case all 
of the non-Gaussianity parameters will have some scale dependence unless the curvaton has exactly 
a quadratic potential, in which case it can be treated as a free test field during inflation. This is 



manifest from eqs. (4.21) and (4.22)"^. In W^ we computed nj^^ for curvaton models with a quartic 



self-interaction term, V — TO^cr^/2 + Acr^, finding a scale dependence proportional to r^o-, which tends 
to be too small to be of observable interest. The result might be different for other type of interactions 
and it would be interesting to compute n/^^ and Ug^^ for generic interacting curvaton models. This, 
however, requires a numerical study and is beyond the scope of the current work |59(| . Here we will 
instead consider the modulated reheating scenario as an example of isocurvature single field models. 
In this case there is little constraint on the form of the modulaton potential and we can use results 
derived in the literature to compute the scale dependencies. 



4.1.2 Modulated reheating 

In this scenario, an isocurvature field a during inflation modulates the decay rate of the inflaton field 
into radiation. Because the expansion rate of the universe changes after the decay, this process can 
convert the initial isocurvature perturbations of the modulaton field into the primordial curvature 
perturbation pi, p7| . This is closely related to the model of modulated preheating pO| and modulated 



trapping |6l|] (see also |62|). This process leads to some level of non-Gaussianity, which depends on 
the efficiency of the transfer, on the functional form of the decay rate r{a) and on the potential of the 
modulaton field V{a). The form of the inflaton potential during horizon crossing is unconstrained, 
assuming the infiaton perturbations can be neglected, but its shape around the minimum does infiuence 
reheating and we assume it has a quadratic potential while it is oscillating. 

For simplicity we will consider the case that F <^ He, where He is the Hubble parameter measured 
at the end of inflation tg . Hence we are assuming that the infiaton decays long after the end of inflation. 
In this case, the curvature perturbation in real space can be written as p7| , H2] 

C(t/,x)^-i^<5a(t„x) + i('-i^] xSa{U,^f + U-^^\ xSa{U,^f + --- , (4.24) 



where ti is a time soon after the horizon crossing of modes of interest. Using Eq. (4.20) we find the 



^For a quadratic model Naaa = (when working to first order in r = Pct/(ZH'^), i.e. considering the curvaton as a 
test field) and h ence (jnt. = 0. As explained in Sec. (pj), in this case we define Ugj^^ = instead of using the formally 



divergent result (4.22) 
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constant parts of /nl and ^nl are given by 

/nl - 5 f 1 - 



rr„ 



5NL 



50 

y 



rr. 



2-3- 



r2r„ 



r^ 



n. 



(4.25) 
(4.26) 



The scale dependencies of /nl and 5nl can be computed using Eqs. (4.21) and (4.22). From these 
equations it is obvious that a potentiaUy large scale dependence, accompanied with large values for 
/nl and ^nl, can arise only if the modulaton field a has large self interactions. For the rest of this 
section we will consider the case of a quartic potential 

A ,, 



V{a) 



4! 



A>0 



(4.27) 



In keeping with the previous literature, we neglect the energy density of the a field after the end of 
reheating, as studying this goes beyond the realms of this project. 



Using Eqs. (|^) and ( [4.22D we find 
5 r Act, 



»^/n 



/nl To 



3H- 



0.1- 



/nlT', 



1/2 
C 



'Snl 



/NL „ 
'SNL 



'Vn 



50 r^ A 



A 



!^n/,, +4x10-3 — 

5nl 5nl /-^c 



(4.28) 



(4.29) 



where rja- = Xaf /{6Hf). In the expression for ng^^ '^c have neglected the contribution proportional 



to r]"^ in Eq. (4.22) which is negligible compared to X/{3Hf) because Xaf <C Hf by construction. 



The quantities nj^^ and 



"■Shl 



could be large, by making a suitable choice of the parameters. 



At first sight, it seems easy to obtain values for these parameters of order 10"^, large enough to 
be detectable, and at the same time compatible with the assumptions that underlie our analysis of 
Section ^ This is correct, but we have to ensure that the parameters satisfy stringent constraints in 
order to obtain acceptable values for the tilt of the spectral index. Indeed, assuming inflation lasted 



considerably longer than 60 efoldings, a natural initial value for the field u is |63 



Ai/4 



(4.30) 



(a different argument changes the power of A from 1/4 to 1/3 and the numerical factors ||6J], but 
the difference is not very important here). Plugging the previous estimate in the expression for 
77ct = Acr|/(67J,j^), and requiring that this parameter is less than 10~^, we find the following bound for 
the coupling A: 



Vcr 



Utt^ 



1/2 



5,10- 



A < 10" 



(4.31) 



The condition F ^ He can place further bounds on 770- since the modulaton is assumed to remain 
nearly frozen until the inflaton decay. We will not further address this issue here. 

Plugging the previous results in (4.28) and (4.29), and using V^ = 2.5 x 10~^ for the normalization 
of the power spectrum, wc find 

600 ^ 20 



P/n 



A^/ 



'Snl I 



I/nlI 
2 X 10^ 

IsnlI 



< 



< 



I/nlI ' 
2 X 10"* 

IsnlI 



(4.32) 
(4.33) 
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where the inequahties are saturated for A ~ 10 . In the estimate for ngp^^j '^^ have neglected the 
first term in Eq. ( 4.29 ), 

9 /nL I I , 3/4 10|/nl| fA -iA^ 

which is subdominant compared to the second term if A^^^ J^ 5 x 10^"'|/nl|- For |/nl| ~ 100, this 
corresponds to A > 6 x 10~^. In the opposite case, A"^/^ ^ 5 x IO^^'I/nlI, the estimate for Ug^^ is given 



by Eq. (4.34) instead of Eq. (4.33) 



We conclude that both nj^^ and Ug^^ could acquire relatively large values, even if the values of 
/nl and (/nl saturate their current observational bounds. It is however important to emphasize that 
I'^/nlI or I'^snlI ^ 0-01 is outside the regime of validity of our formulae, since the accuracy of the 
expansions performed in Sec. starts to become inadequate. 

4.2 Two-field local case 

As a last example, we briefly discuss the so called two field local case, for which f„ and y^ are 



independent of k. This demonstrates an explicit violation of the relation (3.7). In this scenario, the 
formulae at the beginning of this section provide 



6 . V I 



TNL = 7/nl — , (4.35) 

which shows that, in principle, the parameter Wa is an observable. The scale dependences of the 
non-linearity parameters satisfy the following relation 

3 

"tnl = '^SNL = 2'^/nl • (4-36) 

So, as previously stated, we have a different consistency relation between ri/^L and Ut^-^ in this case 



compared to the single field case, Eq. (3.7). Furthermore there is an additional consistency relation 
from n^NL ■ However two-field local models are likely to arise from a test field with a quadratic potential 



1 14 , in which case the amplitude of ^nl tends to be too small to be observable. 

As an explicit example we consider the mixed infiaton-curvaton scenario, assuming the curvaton 

field has a quadratic potential. We discussed this model previously at the level of the bispectrum in 

14 , and found that in a natural limit ti/^l = ^2(n^ ^1)- It therefore follows that for this model rirNL 



is even larger, 

"tn. = -3(nc - 1) . (4.37) 

In this model (/nl ^ /nl [pSj which is too small to be of observational interest |57| . 

5. Shape dependence 

In the previous sections, we concentrated our analysis on the scale dependence of equilateral figures 
(triangles and quadrilatera) . Moreover, we only considered the possibility of varying simultaneously 
all of the sides of the figure by the same proportion. In this section, we study more general situations 
in which scale dependence can arise in parameters characterizing local non-Gaussianity. In particular, 
we consider the case in which the figure under consideration is not equilateral, and the case in which 
we vary the size of only one side, keeping the lengths of the other sides fixed. 
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We start by studying these issues for the parameter /nl, generahzing the arguments developed in 



Sec. 2^, and using the same quantities introduced there. Expanding Eq. ( 2.22 ), around a pivot scale 
fcp using Eq. (2.13), we obtain 



/NL(fcl,fc2,fc3) 



ab 



""3 ( '^multi.o in j.2 



T'f,ab 



Inf^ 



2 perms 



kl 



kl 



kl 



(5.1) 



As in the previous sections, the result is given to first order in \n{ki/k-p) and both the scale-dependent 
and scale-independent parts are given to leading order in slow roll. 



In this approximation, setting for simplicity fcp = 1, equation (5.1) can be re-expressed in a more 
elegant way as 



/NL(fcl,fc2,fc3) 



Y^ rabihk^ 

/ , /nl 

ab 



k 



3+n/,, 



2 perms 



kl 



Z-3 I Z^3 

h.2 -T ^-3 



(5.2) 



Indeed, since both nmuiti.a and rif^ab, for each a, b, are proportional to slow-roll parameters, an expan- 
sion of Eq. (5.2) at first order in slow-roll provides Eq. ( f).l| ). For general single field models, it reduces 
to 



/NL(fcl,fc2,fc3^ 



~ J] 



kV 



Vnl 






Vnl 



+ k. 



S+n 



/nl 



NL 



kl 



k^ 



jL3 



(5.3) 



where /^^ denotes 5fa-a-/S evaluated at the pivot scale and ri/^L = "-/.o-ct- These simple ways of 
expressing the parameter /nl are particularly suitable to analyze how the triangle shape affects the 



scale dependence. The single field expression (5.3) is equivalent to the analogous result given in Section 
3.3 of p4| |, as one can easily check using Appendix M. Eq. (p^ however takes a much simpler form 
than the result in |l4| as a consequence of cancellations that occur when explicitly writing out the 
results in terms of slow roll parameters. 



We note that, although (xE) is not of the form /nl oc (^1/^2^3)"^^'^ which ]1^ used in order 
to make observational forecasts for riy^^: the bispectrum is a sum of three simple, product separable 
terms 



B^{ki, k2, fca) oc {kik2r<-%''''- + 2 



perms , 



(5.4) 



and that it only depends on one new parameter n/„L ■ In the multiple field case the bispectrum will 
typically depend on more parameters than just rif^j^, see (5.2). An exception is the two-field local 



model discussed in Sec. [l.2| , in which case (we also use Eq. ( |4.12| )) 

B^{ki,k2, ks) oc (/fcifc2)"':+("^NL/2)-4 + 2 perms 



(5.5) 



Notice that it therefore follows from ( p.4| ) and (x5) that models with the same /nl and n/^^ can 



have different bispectral shapes which generalise in different ways the local shape. It is possible that 
observations may distinguish between these shapes and that we could therefore learn whether ti/^l 
arises due to single or multi-field effects (or a combination of the two)"*. 

In Sec. 0, we limited our considerations to the scale dependence of /nl for equilateral triangles. 
On the other hand, by means of Eq. (^^) one can observe that, considering a common rescaling for 
all the three vectors, say k^ — > ak^, our previous results remains valid regardless of the triangle shape. 
Namely, 

9 In /nl (afci , afc22_afc3) 



d\n( 



7 ^ /nl (2?^multi,a + nf,ab) 



(5.6) 



CB thanks Sarah Shandera for pointing this out. 
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which is exactly our previous result. 

While the scale dependence of /nl, when simultaneously varying the triangle sides, does not 
depend on the triangle shape, there are other situations in which it does. We might indeed be 
interested on the scale dependence of /nl, when varying the size of only one of the triangle sides, 
keeping the other two fixed (and the triangle closed). In this case, the result does depend on the 
triangle shape. We focus for simplicity on the single-field case, for which the analysis is particularly 
simple, we do not expect our results to change much when considering multiple fields. When varying 
ki — >• aki, equation (5^) becomes 



/NL(afcl,fc2,fc3) 



/nl 



^3fc3 



jL3 

"-2 



l3 



(5.7) 



It is clear that the dependence on a, in this expression, goes to zero in the limit in which fci vanishes. 



This is because the coefficients of the terms depending on a, in equation (5.2), become very suppressed 
with respect to the remaining terms. This situation corresponds to a squeezed triangle: for this shape, 
we then learn that the value of /nl does not change when varying the length of the triangles shortest 
side. 

In order to determine the triangle shape that leads to maximal scale dependence, one is then lead 
to focus on the opposite limit. That is, on configurations for which kf is as large as possible, with 



respect to fc^ + fcf. In this case, indeed, the coefficients of the terms depending on a, in equation ( |3.2[ ), 
become dominant with respect to the other terms. 

This expectation is correct, as shown by the following more detailed analysis. Taking the loga- 
rithmic derivative of /nl along a, we find, at leading order in slow-roll: 



c'ln/NL(a^i,fc2,fc3) 



9 In a 



Vn 



1 



yO 



1- 



3a;'^ Inx + 3y^ \ny 



1 



r 



(5.8) 



where we have defined x = k2/ki, y — k^/ki. We have checked that the terms inside the parenthesis 
are not important for determining the location of the maxima of the previous expression. The maxima 



of Eq. ( p. 81) are therefore determined by the prefactor {1 + x^ + y^)~^ , which is maximized for triangles 
that minimize the combination x^+y^. This corresponds, as anticipated from our previous expectation, 
to the shape for which kf is as large as possible, with respect to the combination fcf + ^1- Calling 9 
the angle between fci and ^2, we have y^ = (1 — x) -I- 2x (1 — cos 9). So we can write 



y 



x^ + (l-x) +2a;(l-cos( 



(5.9) 



It is easy to see that the previous expression is minimized for 9^0 and x ~ 1/2, that is for a folded 
triangle for which k2 — k^ — ki/2. Plugging these values in (^.8|), we find 



a In/ 



NL 



91n( 



1.1 n/„ 



(5.10) 



so we learn that, for the shape that maximizes the scale dependence, we gain around ten per cent 
with respect to the case in which we vary simultaneously all the sides of the triangle. Plots in Fig. |^ 
represent the logarithmic derivative of /nl along a, and graphically show the results discussed so far. 
Notice that the shape which maximizes the scale dependence is indeed given by folded triangles. 



A similar procedure, that generalizes what we have done in Sec. 2.3, can be applied to analyze 
(7nl and tnl- Expanding Eqs. ( ^.29| ) and ( 2.30 ) around a pivot scale fcp using Eqs. ( ^.13 ), ( 2.14 ) and 
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Figure 1: Behavior of the quantity d ln/NL/91na, as a simultaneous function of x (taken between and 
1.5) and of cos 6'. The two plots represent the same figure from two different points of view, that emphasize 
respectively the dependence on x and on cosS. We have chosen rij^^ = 0.01. 



( P.19| ), we obtain the results 
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(5.12) 



We can then proceed with arguments very similar to the ones developed for /nl- Writing k^ — >■ aki 
in Eqs. (5.11) and ( p. 12 ), taking a logarithmic derivative with respect to a and finally setting a = 1, 
we immediately recover the results ( 2.33 ) and ( ^.34 ), derived in Sec. g for n^-^L and Ug^^ for equilateral 
configurations. This shows that for the class of shape preserving variations, k^ — > ak^, the results are 
independent of the figure shape. 

In the single field case, if we vary only one of the sides, say the one labeled by fci, then the scale 
dependence vanishes when fci — >■ 0. We numerically analyzed for which shapes the scale dependence is 
maximal. For the case of f/NL, the analysis is a straightforward generalization of what we did for /nl- 
The shape associated with maximal scale dependence corresponds to a folded polygon, in which three 
of the sides lie over the side whose length is varied. That is. 



fci = 3fc2 = Sfca — 3k. 



A ■ 



(5.13) 



Again, for maximal scale dependence we gain order ten per cent with respect to the case in which we 
vary all the sides simultaneously. 

We also performed a numerical analysis to study tnl , finding again maximal scale dependence for 



the folded shape of Eq. (5.13). For this parameter, we gain around 20 ~ 25 percent with respect to 
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the case in which we vary all the sides by the same amount (this result resonates with the consistency 



relation (3.7)). 



6. Curvature perturbation in coordinate space 

An additional important feature of our approach to the scale dependence of local non-Gaussianity, is 
that it allows one to express the results in coordinate space. In this section, we show how the scale 



dependent coefficients appearing in the momentum space expansion of curvature perturbation ( 2.16 ) 
manifest themselves in coordinate space. It is clear that scale dependence will cause deviations from 
the local form, for which C(x) can be expressed as a power series of a Gaussian variable C^(x) with 
constant coefficients. Using the general single field case as an example, we work out the expression for 
C(x) in the scale-dependent case and quantify how it deviates from the local form. 
In the general single field case, Eq. { p,.9\ ) can be written as 

Ck = Ck + ^/NL(fcp) (^1 + uf^, e{h - k) ln|^ (C^ * C'')k (6.1) 

+ ^5NL(fcp) (l + 7Vl 0{h - k) In- j (C^ * C"^ * C'^)k + • • • , 



which coincides with Eq. (3.1) up to slow roll corrections for constant terms. This form is useful for 
our analysis since the horizon scale ki = aiHi appears explicitly. We have inserted the theta functions 
9{ki — k) to explicitly indicate that the result holds only for super- horizon modes k < hi. ki > kp 
should correspond to a physically smaller scale than any of the modes of interest. Recall that similar 
theta functions are included in our definition of (^, Eq. (2.10|). Therefore Ck' can be viewed as a 



smoothed quantity; in Fourier space the window function is simply a top hat with the cutoff set at 
the horizon scale fc,. 



Taking the inverse Fourier transform of (6.1) we find 

^/NL(fcp)C«(x)^-H| 



C(x) = C''(x) + ^/NL(fcp)C°(x)2 + _L5NL(fcp)C''(x)3 (6.2) 



3 /■ dk k 

+ 5/NL(fcp)n/,, J ^^e*-0(fc, _ fc)(C° ^C^)kln- 

q /■ dk k 

+ ^9Mkp)n,^.J ^^e*-0(fc,-fc)(C«K^K%ln- + ... . 

The two integrals describe deviations from the local form. They can be written more explicitly by 
performing the following manipulations 

"^^ e^k.x^(fc,-fc)(C°^C^)kln^ = /dy f -^e^M--y)0ik, ^ k) (""{yrin^ (6.3) 



(27r)3 ' ' "^ ^ "^ h J J (2^)3 ^ ' '^ '■" k 

/"j^^G/ ^2 1 sin(fc,|x-y|)-Si(fc,|x-y|) 

= /dyC°(y)^/(|x-y|), 
and similarly for the second integral. Using this we can rewrite Eq. (p.2|) as 



C(x) = C^(x) + -/NL(fcp) ( C^(x)^ + nj,, I dy/(|x - y|)C^(y)^ ) (6.4) 

+ 4ffNL(A:p) (C^(x)3+n,,, / dy/(|x-y|)C^(y)3 ) + 



19 



This result clearly shows how the scale dependence of faa and g^aa in Eq. ( ^.16| ) renders C(x) a 
nonlocal function of C*^(^)- Because of the integrals in Eq. (6.4), the curvature perturbation C(x) can 



not be expressed in terms of C^'^ (x) evaluated at the same point x but one needs to know C,^ (y) in the 
entire region where /(|x — y|) is non- vanishing. The behavior of /(|x — y|) is depicted in Fig. o which 
also displays the inverse Fourier transform, 



m|x-y|) = 



sin(fc^ |x - y|) - fc^ |x - y| cos(fc^ |x - y|) 
27r2|x-yP 



(6.5) 



of the top hat window function 9(ki — k) included in the definition of Q^, Eq. (2.10). Both W{x) 



|W(x 





Figure 2: Absolute values of the functions W and / plotted on logarithmic scales for the choice fci = 4 (in 
arbitrary units). 



and I{x) are approximatively constant at scales x < [ki)^^ . They both fall off for x > (ki)^^ but 
I{x) remains negative definite unlike W{x) which starts to oscillate rapidly. Keeping in mind that a 
smoothing over W{x) is implicit in the definition of ('~^{x), we therefore see that the convolutions of 
^ (x) with I{x) in Eq. (p^) pick up a non-trivial contribution from the superhorizon modes x > {ki)~^ 
where W{x) effectively falls off faster than I{x). This contribution makes Eq. ( |6.4[ ) deviate from the 
local form. 

The analysis can in principle be straightforwardly generalized to the multi- field case. The difference 
compared to the general single field case is the appearance of several unrelated Gaussian fields C] 



G,a 



in Eq. (2.16). This in general makes it impossible to write C(x) as a series of a single Gaussian field 
even if the coefficients in Ec^. ( [2.16 ) would be constants. 



7. Conclusions 

We have discussed a new approach, based on the (5iV-formalism, for studying the scale dependence of 
non-Gaussianity parameters. We have obtained explicit expressions for the scale dependence of the 
quantities /nlj ''nl E^nd i^nl associated with the bispectrum and trispectrum of primordial curvature 
perturbations. Our results depend on the slow-roll parameters evaluated at horizon exit, and on the 
derivatives of the number of e-foldings and the inflationary potential. The parameters controlling the 
scale dependence of non-Gaussianity depend on properties of the the inflationary potential, namely its 
third and fourth derivatives, which in all observationally interesting cases cannot be probed by only 
studying the spectral index of the power spectrum and its running. 

As a consequence, the scale dependence of non-Gaussianity provides additional powerful observ- 
ables, able to offer novel information about the mechanism which generates the curvature perturba- 
tions. We demonstrated these features in the concrete example of modulated reheating. In models 
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with a quartic potential for the modulating field, we have shown that the associated non-linearity 
parameters, and their scale dependence, can be large enough to be observable. 

While in most of the discussion we worked in momentum space, in the last part we also discussed 
how to describe our results in coordinate space. We provided an expression for curvature perturbations 
in coordinate space, that generalizes the frequently used local Ansatz, and that exhibits directly in 
real space the effects of scale dependence of non-Gaussian parameters. 

Our results allow us to put onto a firm basis the phenomenological parameterizations of the scale 
dependence of non-Gaussian observables. In many models of observational interest, our formulae are 
relatively simple and depend on a single new parameter, the scale dependence of the non-linearity 
parameter. It would be interesting to use these results for analysing or simulating non-Gaussian 
data. At the same time, our investigation allows us to identify which properties inflationary models 
have to satisfy, in order to obtain large non-Gaussianity with sizeable scale dependence. It would be 
interesting to apply it to analyse further models, for example those in which multiple fields interact 
during infiation or where the non-Gaussianity is generated by an inhomogeneous end of inflation. 
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Appendices 

A. Explicit expressions for ri/at and rig at 

For a system of slowly rolling scalar fields (j)a, the equations of motion are given by iHcjia — ~Va and 
3H^ = V, to leading order in slow roll. Here we are interested in the evolution during a short time 
interval from to to ti. The slow roll equations can easily be solved for (/)a as 

nU) = rih) - V2^\n{aJao) + O (^e^/^ln'iajao)) , (A.l) 

where we have used the identity Hdt = din a, which holds to leading order in slow roll. In the following 
we use the notation ©(e") to denote the combinations of the slow roll parameters Eq, rjab of order e". 



Differentiating Eq. (A.l) with respect to 4'°' (to) and keeping the number of e- foldings ln(ai/ao) 



fixed, we can compute the coefficients appearing in Eq. (2.4). We choose the initial time to as the time 
tk of horizon crossing of a mode fc, defined by a^Hk — k. Using hi{ai/ak) = ln{aiHi/k), which is valid 
at leading order in slow roll, the three first coefficients in Eq. (^^) can be written as 



^rit^) 



= Sa 



dcji^tk) 



d4\tk)d(j,-{tk) 

d^'ritj) 

dcp>^{tk)dr{tu)dcp'i{tt,) 




F\' In 

aoc 



F^^^ In 
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where the primes in 0{V"') etc. denote derivatives with respect any of the fields (f>a, and 



abc 

r^(3) 
abed 



V2 i -i^eaEbec + r]ab\^ + r]bcV^ + VcaV^ - 



1 V, 



abc \ 



^/23iJ2; ' 



--iiVadV^b^ + ribd\/<^cf^a + VcdV^aeb) 



(A.5) 
(A.6) 



-^/e^ ( 24^eaebec - 'irjabV^ - '^r\bc\f^a - '^rica^f^ + ^2^-^ 



I2l. 



■K, 



bd 



1^. 



-Vbcd 



.3^2 ' V" a3^2 

Substituting these into Eq. ( 2.f5| ), 



/2^ 



^^« 



3il2 



3i?2 

'7ab'7crf + ??fcc?7ad + VcaVbd 



Vn 



bed 



3iJ2 



"^g.abc 






iVa, 



f)C 



d6c 



^d ^(3) 



A^a. 



fcc 



dabc 



(A.7) 
(A.8) 



we obtain fully explicit results for the parameters nf,ab and Ug^bc- The results are derived retaining 



only terms up to first order in hi{aiHi/k) in Eqs. (A. 2) - (A. 4). Higher order terms are suppressed 
by slow-roll parameters and their combinations with the derivatives of the potential. The former are 
small by construction and the latter also naturally remain small, provided that the flatness of the 
scalar field potential during infiation is not a result of extreme fine-tuning. Furthermore, since the 
logarithms never grow very large for the observable super-horizon modes, \n{aiHi/k) < 0(10), we 
conclude that the higher order contributions can indeed be neglected at first order. 

B. On the different formulations of the 5N approach 

The 5N expression for the super-horizon-scale curvature perturbation 



dq 



Ck(t/) - Y. Naitf,u)srwiu) + ^Y. ^-b{tf, u) f jp 

^ ab •' ^^'^ 



iq(iO'^4-q(i.) + 



(B.l) 



is by construction independent of the choice of the initial spatially flat hypersurface ti > tk- This 
property follows from the definition of the curvature perturbation |Q as demonstrated in iQ. A 
commonly used choice is to set ti equal to the time tk of horizon crossing of the mode k. The analysis 
in [|l4| was performed using this choice. In our current work, we have instead chosen ti as a time soon 
after tk following e.g. |23]. Here we explicitly compare the two choices. For simplicity, we consider 



only terms up to second order in Eq. (B.l). Generalization to higher orders is straightforward. 



Using the chain rule it is easy to switch between the coefficients Nab.Xifi^i) ^^d Nab..itf,tk) in 
the two different formulations. For the first and second order terms shown in Eq. (|B.l|) we obtain 



Naitf,t, 
Nab{tf,ti 



^ dN{tf,t,) 
driu) 



E 



dN{tf,tk)d(j)''{tk) 



b'itk) driu) 



d^N{tf,t,) 



E 



dN{tf,tk) d^^itk) 



dr{U)dcf>b{t,) ^ d(^-(tk) d(l)-{ti)d4>b{ti) 
ntk)d<P''{tk) d^N{tf,tk) 



(B.2) 
(B.3) 



E 

cd 



driu) dnti) dritk)dntk) 
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In the second equality of both equations we have replaced the time argument U in N{tf, ti) by tk making 
use of the fact that ti and tk both label spatially flat hypersurfaces. This implies that N{ti,tk), the 
number of e- foldings from tk to ti, is a constant under differentiation with respect to the fields. Writing 
N{tf,ti) = N{tf,tk) — N{ti,tk), we thus immediately see that N{tf,ti) can be replaced by N{tf,tk) 
in Eqs. (^) and (j^J). 

Substituting Eqs. ([B.2|) and (B.3) into Eq. (B.l), we obtain 



Ck(t/) = ^A^a(t/,tfc) 



d^itk) 



1 d'ritk) 



b{t,) "-^'^^^'^ ■ 2 dcp'>{t,)dnt^) 



— --350^(ii)(50S^_q(i,) 



(B.4) 



lY.^-b{tf,tk) 



dritk) d4>\tk) 



''{U 



\U 



-^,mi^¥'i^i-.{t^) 



On the other hand, according to Eq. (2^) we have 



Hlitk) 



d(i>\t,) ' 



t>iiu 



dq 



1 d^rjtk) 

2d(j)\t,)dc^-{t,)J (27r 



■54>\{U)5c^l^^{t^) 



(B.5) 



In arriving at this result we have first taken the Fourier transform of Eq. (2^) and only thereafter set 



one of the time arguments equal to tk- Using Eq. (B.5) we can rewrite Eq. (B.4) as 



a ab J \ ) 



k)S4'i-jtk) 



(B.6) 



This way of writing Ck;(i/) is equivalent to Eq. (BT) and the two expressions differ formally only by the 
choice of the initial time ti , as expected. The relation between the coefficients in the two formulations 
is given by Eqs. ( [B.2| ) and (B^), and the field perturbations are related by Eq. (B_^). These results 
explicitly show how to switch from one formulation to another. 

In iQ, the result for Uf^^, measuring the scale dependence of /nl, was expressed in terms of 
the parameters nj = din Ni{tf ,tk) /d\n k and nij = d\nNjj{tf,tk)/dlnk, see e.g. Eq. (69) in that 
paper. (Here we follow the notation of (|lj| and label the scalar field species ^/ by capital letters. This 
also serves to distinguish the parameters nj and njj from the quantities defined in the current work.) 
Using Eqs. ( |B.2| ) and (B_^) together with the results derived in Appendix Q we readily obtain explicit 
expressions for these parameters 



ni = 



nij 



d\n Ni{tf,tk) _ 
d In fc 

d\nNij{tf,tk) 
dlnfc 






rifJJ 



Ef% 



K 



V Nij 



Nij 



(B.7) 
(B.8) 



In the rightmost expressions we have suppressed the time arguments i,; for brevity, e.g. Nj = Ni(tf, ti). 
Using these results, it is straightforward to check that the general expression given for rif^^ in Eq. (69) 
of agrees with our Eq. (|2.24D. 
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